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The energy dependence of the observable two-particle correlator in 
search for local strong parity violation in Au+Au collisions is estimated 
within a simple phenomenological model. The model reproduces available 
RHIC data but predicts that at LHC the chiral magnetic eifect (CME) will 
be about 20 times weaker than at RHIC, contrary to the first LHC mea- 
surements. In the lower energy range this effect should vanish sharply at an 
energy somewhere above the top SPS one in agreement with the preliminary 
results of the Beam Energy Scan program. To elucidate CME background 
effects a transport HSD model including magnetic field evolution is put 
forward and electromagnetic dynamics at RHIC energies is investigated. It 
is observed that the electromagnetic field included into the hadronic model 
does not influence on observables due to mutual compensation of effects of 
electric and magnetic fields. 

PACS numbers: 25.75.-q, 25.75.Ag 

1. Introduction 

The existence of nontrivial topological configurations in QCD vacuum 
is a fundamental property of the gauge theory though until now there is 
no direct experimental evidence for topological effects. Transitions between 
different topological states occur with the change of the topological number 

characterizing these states and induce anomalous processes like local 
violation of the V and CV symmetry. The interplay of these topological 
configurations with (chiral) quarks results in the local imbalance of chiral- 
ity. Such chiral asymmetry coupled to a strong magnetic field, created by 
colliding nuclei perpendicularly to the reaction plane, induces a current of 
electric charge along the direction of a magnetic field which leads to a sep- 
aration of oppositely charged particles with respect to the reaction plane. 
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Thus, as was argued in Refs. [H El O [U [5] the topological effects in QCD may 
be observed in heavy ion collisions directly in the presence of very intense 
external magnetic fields due to the "Chiral Magnetic Effect" (CME) as a 
manifestation of spontaneous violation of the CV symmetry. It was shown 
that the electromagnetic field of the required strength can be created in 
relativistic heavy-ion collisions O [6] . First experimental evidence for the 
CME identified via the observed charge separation effect with respect to 
the reaction plane was presented by the STAR Collaboration at RHIC [7j. 
Recently preliminary data below the nominal RHIC energy (Beam Energy 
Scan program) [8j and at the LHC energy [9] have been obtained. 

In this paper we shortly remain our analysis of the first STAR data jlOj 
and compare results with new finding. In respect to the CME, we also 
make a dynamical estimate of the CME background based on the nonequal- 
ibrium Hadron-String-Dynamics (HSD) microscopical transport approach 
[llj supplemented by the treatment of electromagnetic field evolution. 



We remind here our phenomenological model used in the CME analysis 
in Ref. [10] . 

For one-dimensional random walk in the topological number space the 
topological charge (winding number) generated during the time tb, when 
the magnetic field is present, may be estimated as 



where F^ is the sphaleron transition rate which in weak and strong coupling 
Ts ~ with different coefficients. The initial temperature Tq of the pro- 
duced matter at time r ~ l/Qs is proportional to the saturation momentum 
Qs, Tq = c Qs- In the last semiequality of ([T]) the expansion time and the 
corresponding time dependence of the temperature are neglected. Since siz- 
able sphaleron transitions occur only in the deconfined phase, the time tb 
in Eq. ([1]) is really the smallest lifetime between the strong magnetic field 
Tb one and the lifetime of deconfined matter Tgi tb = minify, r^}. 

The measured electric charge particle asymmetry is associated with the 
averaged correlator a by the following relation |12j : 

(cOs(V'a + 1pl3- 'i'^Rp)) = (C0s(^a + # " 2^c))/'f2,c = Vi^aVl,l3 - a^ajB ,(2) 

where ^ rp is the azimuthal angle of the reaction plane defined by the 
beam axis and the line joining the centers of colliding nuclei. Averaging 
in is carried out over the whole event ensemble. The second equality 
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in ([2]) corresponds to azimuthal measurements with respect to particle of 
type c extracted from three-body correlation analysis [12], vi and V2 are the 
directed and elliptic flow parameters, respectively. According to Ref. [Ij an 
average correlator a = ^aaap is related to the topological charge, 

n,„ pTB I r- ^-1/16 r— fo\ 

JT- ~ W^NN) ' ■ VTB, (oj 



f^A^hadrons/'^y VdA^hadrons/d?/ V 



where absorption and rescattering in dense matter are neglected for the 
same and opposite charge correlations. In the last equality we assumed 
that ~ s]^^! ~ dN^sdrons/dy [13J. 

As follows from Eq. Q, the CV violation effect can be quantified by the 
correlator 

= Kau iV^y^^^ ■ . (4) 

The normalization constant K^u can be tuned at the reference energy 
^i]^ =200 GeV from the inverse relation (|3]) and experimental value dexp 
at this energy for b =10 fm 



a 



2 



• (200)1/8 



" "rB(200) • 

The only quantity needed is the lifetime tb defined as the time during 
which the magnetic field is above the critical value Bcrit needed to support 
a fermion Landau level on the domain wall eBcrit = '^'^/Sd, where is 
the domain wall area. Since the size of the domain wall is not reliably 
known, it is hard to pin down the number, but it should be of the order 
of m^. Thus, we have to treat Bcrit as a free parameter which defines 
the lifetime from the calculated B{t) distribution [TU]. As to the impact 
parameter distribution, the CME is assumed to be roughly linear in b/R 
Ref. f3]. Taking this as a hypothesis we evaluate the centrality dependence 
of the CME fitting this hue to points 6 = 10 fm (or centrality (40 - 50)%) 
to be estimated in our model and 6 = where the CME is zero. The results 
are presented in Fig.[T]for Au+Au collisions. As is seen, the calculated lines 
quite reasonably reproduce the measured points of azimuthal asymmetry of 
same-charge particles for Au+Au collisions at ^snn =200 and 62 GeV. The 
chosen critical field value eBcrit = 0.7 results in absence of the CME 
above the top SPS energy because the critical magnetic field practically 
coincides with the maximal field at this bombarding energy. This finding 
is in a qualitative agreement with very recent preliminary STAR results 
that the difference between same-charge and opposite charge correlations 
is decreasing with decreasing beam energy what takes place at ^/snn <40 
GeV [S]. 
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Fig. 1. Centrality dependence of the CME. Experimental points for Au+Au colli- 
sions are from [7]. The dotted line is our prediction for Au+Au collisions at the 
5.5 TeV energy. 

In the model considered, the CME at the energy 5.5 TeV is expected 
to be less by a factor of about 20 as compared to those at the RHIC en- 
ergy [8]. Note that at the LHC energy we applied a simplified semi-analytical 
model [6] for magnetic field creation. Thus, our LHC estimate should be 
considered as an upper limit for the CME. Recent preliminary LHC PbPb 
data [S] show a remarkable agreement in both the magnitude and the behav- 
ior with the results reported by STAR in Au-Au collisions at ^/SNN =200 
GeV. Preliminary kinetic calculations at the LHC energy 2.76 TeV increase 
the maximal magnetic field by a factor about three but it does not influence 
essentially on the relaxation time. 



The discussed CME signal, the electric charge asymmetry with respect 
to the reaction plane, can originate not only from the spontaneous local CP 
violation but also be simulated by other possible effects. In this respect it is 
important to consider the CME background. We shall do that considering 
a full evolution of nucleus-nucleus collisions in terms of the HSD transport 
model [TI] but including formation of electromagnetic field as well as its 
evolution and impact on particle propagation. 

Generalized on-shell transport equations for N strongly interacting par- 
ticles in the presence of electromagnetic fields can be written as |l4j 



-+(eE + {-)vxB)VpA f{r,p,t) = Icoii{fJi,:.fN-i), (6) 
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which are supplemented by the wave equation for the magnetic field whose 
solution in the semi-classical approximation for point-like moving charges is 
reduced to the retarded Lienard-Wiechert potential [6]. The quasiparticle 
propagation in the electromagnetic field is calculated according to dp/dt = 
eE + (e/c) V X B. 

In a nuclear collision, the magnetic field will be a superposition of 
solenoidal fields from different moving charges. It is illustrated in Fig. [21 
Its first panel is taken at a quite early compression stage with t =0.05 fm/c 
close to the maximal overlapping where the magnetic field here is maximal. 
The overlapping strongly interacting participant region has an "almond" - 
like shape. The nuclear region outside this almond corresponds to spectator 
matter which is the dominant source of the electromagnetic field at the very 
beginning of the nuclear collision. Note that in the HSD code the particles 
are subdivided into target and projectile spectators and participants not 
geometrically but dynamically: spectators are nucleons which suffered yet 
no collision. The next time moment corresponding to expansion stage is 
illustrated in Fig. [5] and in more detail studied in Ref. [Hj. The time 



AuAu, = 200 GeV, b=10 fm, t=0.05 fm/c AuAu, V5„„ = 200 GeV, b=10 fm. t=0.2 fm/c 




Fig. 2. Distribution of the magnetic field strength eBy in the y = plane at 
t =0.05 and 0.2 fm/c for Au+Au collisions at ^/snn =200 and b =10 fm. The 
collision geometry is projected onx—z plane by points corresponding to a particular 
spectator position. Curves (and their projections) are levels of the constant eBy. 

evolution of the magnetic field and energy density at the center of the al- 
mond region eBy{x = 0,y = 0,z) for Au+Au collisions for the colliding 
energy y^s^vjv =200 GeV and the impact parameter 6 =10 fm is shown in 
Fig. [3l It is seen that the largest values of eBy ~ 5m^ are reached for a 
very short time. Note that this is an extremely high magnetic field, since 
~ 10^^ gauss. Then, the system expands and the magnetic field de- 
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creases. It is of interest to note that in our transport model, the spectator 
contribution to the magnetic field is practically vanishing at the center for 
t ~1 fm/c. In subsequent times the magnetic field eBy is formed essen- 
tially due to produced participants with roughly equal number of negative 
and positive charges which approximately compensate each other ^14j. The 
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Fig. 3. Time evolution of the magnetic field strength eBy (left) and energy density 
(right) at the central point x =0, y = and z = fm/c for Au+Au collisions at 
various RHIC energies. 



evolution of the energy density of created particles is presented in the right 
panel of Fig. [3l Here the maximal energy density (in the center of the col- 
liding system) is e > 50 GeV//m^ at the moment of maximal overlap of the 
nuclei. When the system expands, it takes a sausage-like shape and the en- 
ergy density drops fast. But even at the time t ~ 0.5 fm/c the local energy 
density is seen to be above an effective threshold of a quark-gluon phase 
transition e >1 GeV//m^. As shown in [14] the location of the maximum 
energy density correlates with that for the magnetic field. 

The background electric field, being orthogonal to the magnetic one, 
is directed along the x axis. The evolution of the eEx field for peripheral 
(6 =10 fm) collisions of Au+Au at the top RHIC energy is presented in 
Fig. m Similar to the case of the magnetic field, the eEx{x,y = 0,z) dis- 
tribution is also inhomogeneous and closely correlates with geometry while 
the field strength looks "hedgehog" shaped. When the two nuclei collide, 
the electric fields in the overlap region significantly compensate each other, 
and the electric field E in the target and projectile spectator parts have 
opposite signs. As a result, the locations of the maximum/minimum are 
not in the central point of the overlap region - as they are for the magnetic 
field - but shifted slightly outside. The maximum of the electric field can be 
quite large. All these features are seen explicitly in Fig. U] where the tempo- 
ral evolutions of eEx{x,0, z) and eEy{x,0, z) are given. Due to destructive 
interference or the "hedgehog" effect, the electric field in the central part 
of the overlap region (x ~ fm) is consistent with zero apart from a short 
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AuAu. \„ = 200GeV, b=10fm, t=O.OI fm/c AuAu, V%=200GeV, b=10fm, t=0.05fin/c 




Fig. 4. (Color online) Evolution of the x- and y-components of the electric field at 
incoming and maximal overlap in Au+Au(-\/s= 200 GeV) collisions at the impact 
parameter b =10 fm. The eE^ =const levels and spectator points are shown in the 
projection on the {x — z) plane. 

period just before reaching maximal overlap. Note that the electric field at 
the central point is negligible for t > 0.15 fm/c. 

4. Observables and electric charge separation 

The HSD model quite successfully describes many observables in a large 
range of the collision energy. Here we investigate to what extent the elec- 
tromagnetic field - incorporated in the HSD approach - will affect some 
observables. We shall limit ourselves to Au+Au collisions at ^snn = 200 
GeV and impact parameter 6 = 10 fm. Here we calculate the whole nuclear 
interaction including the decays of resonances at least up to times of 50 
fm/c. 

The HSD results for the versions without and with electromagnetic field 
are presented in Fig. [5j With a high degree of accuracy, we see no difference 
between these two versions in the transverse mass mt and rapidity y. In 
Fig. [6] the transverse momentum dependence of the elliptic flow of charged 
pions is compared for two versions (with and without field) of the HSD 
model. We do not observe any significant difference between the two cases. 
Slight differences are seen in the range of ~ 1 GeV/c but certainly it 
can not be considered as significant. Note that generally the HSD model 
underestimates the elliptic flow, but an inclusion of partonic degrees of 
freedom within the PHSD approach allows to describe perfectly well the 
Pt dependence of V2 at the top RHIC energy [15j. The correlator ([2]) is 
calculated on the event-by-event basis. The experimental data from the 
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Fig. 5. Transverse mass and rapidity distributions of charged pions produced in 
Au+Au {y/SNN = 200 GeV) collisions at 6 = 10 fm. The results calculated with 
and without electromagnetic field are plotted by the dotted and solid lines, resp. 
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Fig. 6. The transverse momentum dependence of the elliptic flow for Au+Au 
{y/sNN —200 GeV) coUisions at b =10 fm. 



STAR Collaboration [7] and the results of HSD calculations are presented 
in the left panel of Fig. [71 The experimental acceptance \ri\ < 1 and 0.15 < 
Pt < 2 GeV were also incorporated in theoretical calculations. Note that 
the theoretical reaction plane is fixed exactly by the initial conditions and 
therefore is not defined by a correlation with a third charged particle as in 
the experiment 0. As is seen, the HSD model shows no charge separation 
effect. The reason of that is explained in the right panel of Fig. [71 The 
average momentum increment Ap due to electric and magnetic fields is 
almost completely compensated in every component. 

5. Conclusions 

The model energy dependence of electric like-charge pairs can be recon- 
ciled with experiment [7] by a detailed treatment of the lifetime taking into 
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Fig. 7. Left: Azimuthal correlation in the transverse plane vs centrality for like and 
unlike charged pions from Au+Au {yfSNN = 200 GeV) collisions. The experimental 
points - connected by lines - are taken from [7]. Right: Time dependence of the 
momentum increment of mesons created in Au+Au {^/s = 200 GeV) collisions with 
the impact parameter 6 = 10 fm within HSD. Numbers along the right axis show 
the total transverse momentum increment Ap due to electric and magnetic fields. 



account both the time of being in a strong magnetic field and time evolu- 
tion of the energy density in the QGP phase. For the chosen parameters 
we are able to describe RHIC data for Au+Au collisions on electric charge 
separation at two available energies. We predict that the effect will be much 
smaller at the LHC energy and will sharply disappear near the top energy 
of SPS. Coming experiments at the Large Hadron Collider [9] and that of 
the planned Beam Energy Scan program at RHIC [8] are of great interest 
since they will allow one to test the CME scenario and to infer the critical 
magnetic field eBcrit governing the spontaneous local CP violation. 

We have extended the hadron string dynamics model for describing the 
formation of the retarded electromagnetic field, its evolution during a nu- 
clear collision and the effect of this field on particle propagation. The case 
of the Au+Au collision at ^/sN^ = 200 GeV for 6 = 10 fm is considered 
in great detail. It is shown that the most intensive magnetic field oriented 
perpendicularly to the reaction plane is formed during the time when the 
Lorentz-contracted nuclei are passing through each other, t < 0.2 fm/c. 
The maximal strength of the magnetic field here attains very high values. 
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eBy/rri^ ~ 5. However, due to the compensation effect, tfie electromagnetic 
field does not influence on observables and, in particular, on the asymmetry 
of charged mesons with respect to the reaction plane. 
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